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Thermodynamic and spectral properties of a photon gas 
in e'''e~7 plasmas are studied. The effect of a finite effective 
mass of a photon, associated with the plasma frequency cutoff, 
is self-consistently included. In the ultra-relativistic plasma, 
the photon spectrum turns out to be universal with the tem- 
perature normalized plasma frequency cutoff being a funda- 
mental constant independent of plasma parameters. Such a 
universality does not hold in the non-relativistic plasma. 
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It is well known that a photon gas which is in thernio- 
dynamic equilibrium with matter exhibits the Planck- 
ian blackbody energy spectrum. This is correct only 
if matter is electrically neutral (for example, an atomic 
gas) where the thermodynamic equilibrium is established 
via photon emission and absorption in atomic transi- 
tions. Such an equilibrium no longer forms in ionized 
gases (plasmas) where the number of charged particles 
is not negligible. An electromagnetic wave (photon) is 
strongly coupled to collective plasma excitations via os- 
cillating electric and magnetic fields. The dispersion re- 
lation of photons in a plasma is cj^ ~ cjp -I- fc^c^, where 
Wp is referred to as the plasma frequency which depends 
on plasma parameters (e.g., densities of electrons and 
positrons) . Quite interestingly, such a dispersion may be 
interpreted as a manifestation of a finite effective mass 
of a photon, m-y ~ fiwp/c^. There are no photons with 
energies lower than its effective rest energy. Hence, one 
may expect a cutoff in the Planckian spectrum for fre- 
quencies oj < Up. A new, modified blackbody, distri- 
bution of photons forms in this case. The shape of the 
spectrum is controlled by the densities of electrons and 
positrons. It is remarkable that in the relativistic plasma, 
the electron and positron number densities are, in turn, 
self-consistently determined from the equilibrium between 
the photons and the electrons and positrons with respect 
to the e"'"e~-pair production. 

Objects with hot, positron-electron-photon (e+e^7) 
plasmas, referred to as fireballs, are not rare in astro- 
physics. For example, such fireball conditions are be- 
lieved to be realized at early stages of the formation of 
the Universe (see, for example, |^), in energetic cosmo- 
logical explosions which produce 7-ray bursts (see H for 
review), etc.. Hence, the study of thermodynamic and 
statistical properties of a gas of photons in the e^e^7 
fireballs is of both academic and practical interest. 

In this paper, we find a self-consistent, blackbody ra- 
diation spectrum and thermodynamic properties of the 
photon gas in the e^e~7 fireball plasma. In particu- 



lar, the equilibrium spectrum of photons in the ultra- 
relativistic case (T 3> m^c^) turns out to be universal: 

1. the frequency cutoff normalized by temperature is 
a fundamental constant independent of any plasma 
parameters; 

2. the simple, power-law scalings of thermodynamic 
parameters such as pressure, density, etc. vs. tem- 
perature, valid for the Planck distribution, hold for 
"massive" photons as well. Numerical pre-factors, 
however, modify self-consistently. 

In the end of the paper, we also briefly consider a non- 
relativistic case and discuss implications for the cosmic 
microwave background spectrum distortion, purely radia- 
tive fireball model of 7-ray bursts , and stellar interiors. 

Here we should comment on one important point. The 
gas of photons obeys the Bose-Einstein statistics with 
a vanishing chemical potential and, also, is assumed to 
be in thermodynamic equilibrium with matter which, in 
general, has a non-vanishing chemical potential. On the 
other hand, one of the conditions of thermodynamic equi- 
librium is the equality of the chemical potentials of the 
interacting systems. This contradiction may be resolved 
by treating the process of radiation as a first-order phase 
transition (Bose condensation) and properly taking into 
account the equality of the chemical potentials of mat- 
ter and radiation. An excellent discussion of this subject 
may be found in Rcf. [0 . The deviation of a photon spec- 
trum from the Planck distribution may be significant in 
the non-relativistic plasmas, where thermodynamic equi- 
librium is established via, for example, Compton scat- 
tering of a photon on e+ or e~ . In the ultra-relativistic 
regime (considered in this paper) two- and three- annihi- 
lation/creation of an e"'"e~-pair dominate and are respon- 
sible for establishing thermodynamic equilibrium. The 
chemical potential of an e~'"e~-pair is identically zero (by 
definition, e+ and e^ have equal but opposite in sign 
chemical potentials) and the condition of the equality of 
the chemical potentials of matter and radiation is satis- 
fied automatically. The Planck statistics of photons is 
thus valid in the ultra-relativistic e^e~^ plasma. 

Ultra-relativistic case — The dispersion relation of 
electromagnetic waves in the ultra-relativistic plasma 
reads [01: 
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Yla ^T^c^UaC^ I2>T is the relativistic 
plasma frequency, Ca and n^ are the charge and num- 
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ber density of species a, the sum, ^a, goes over all rela- 
tivistic (T ^ niaC^) species, i.e., electrons and positrons 
in our case, and the contribution from non-relativistic 
species being mc^/T times smaller is neglected here. An 
exact analytical solution to this equation is unknown. 
The asymptotic behavior 
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suggests, however, the following approximate solution to 
Eq. (p[D: 
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which roughly satisfies both asymptotics and, also, is the 
exact dispersion relation in the non-relativistic plasma. A 
direct numerical comparison of Eqs. p) and (Uj) shows 
that the fractional error of the approximation does not 
exceed 3%. 

In the ultra-relativistic plasma, the electrons and 
positrons are in thermal equilibrium. Hence, their num- 
ber densities are equal and ne+ = Uf,- — f«7, where the 
factor 7/8 accounts for the difference in Fermi and Bose 
statistics pi and n^ is the number density of photons. 
The relativistic plasma frequency becomes now 
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Note here that the plasma frequency depends on the den- 
sity of photons to be determined |^ . 

The occupation numbers for a gas of photons are given 
by the Bose-Einstein distribution with a vanishing chem- 
ical potential: n(fc) = [exp (?kj(fc)/T) — 1] . The num- 
ber density of photons is formally written as follows [p| : 

T^3 



T^^h^c^ Jo exp (y'lTre'^h'^c-^n^/ST^ + xA - 1 



(4) 

where we made use of the fact that the degeneracy g = 2 
for photons (two independent polarizations). This equa- 
tion self- consistently determines the number density of 
photons in the e+e~7 fireball plasma where the plasma 
frequency cutoff (i.e., the effective photon mass), in turn, 
depends on n^. 

To proceed further we note a remarkable fact. Namely, 
n^ and T enter Eq. (jj) in combination n^/T'^. Then, 
the solution to this equation must be sought in the form 
n^ — const. X T^ with the const, to be determined from 
Eq. (0). Thus, the parameterization n^ ex T^ removes ex- 
plicit temperature dependence of the integral in Eq. (|j)! 
Hence this solution for n^ is universal. Note, this is usu- 
ally not the case for particles with nonzero mass. Since 
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enters the integral over particle momenta, p, via the 
combination {mc^ /T^, so that all the thermodynamic 
quantities (e.g., density, energy, entropy, etc.) do not 
exhibit universal, power-law scalings. For photons, how- 
ever, m^ ex T and the universal (Planckian) scalings are 
recovered. Numerical coefhcients, however, modify self- 
consistently. 



Now, we write n^ 



f npi, where npi 



[2C(3)/7r2] {T/hcf and C(3) ~ 1.2 is the Riemann zeta 
function of 3 and / is a constant to be calculated. Here- 
after the subscript "PI" denotes values calculated for the 
Planckian distribution. Then Eq. (H) becomes 
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and ae — e^ /he is the fine structure constant. Given /, 
the energy density, pressure of the gas of photons as well 
as other thermodynamic parameters may be calculated 
straightforwardly B: 
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Numerically solving Eq. (0) for / = (1 + Sfn) and using 
Eqs. (0), one finally obtains: 

(8a) 
(8b) 
(8c) 

where epi = Spn = (tt^ /15){T'^/h^c^). It is well known 
that the equation of state of a gas of non-relativistic 
massive particles, T <C mc"^, is Pnr = fenr and that 
of ultra-relativistic (or massless) particles, T ^ mc^, is 
Pr — -iei-. Since photons in the relativistic plasma turn 
out to be "sub-relativistic" , fkVp^rei/T ~ 5 x 10"^, we ex- 
pect that the effective equation of state of photons will 
be p^ = ae-y, ^ < a < |. Indeed, from Eqs. (H) it follows 
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The blackbody equilibrium radiation spectrum from 
the ultra-relativistic e+e~7 fireball plasma is easily ob- 



tained from Eq. 

de^(cj,) : 
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where w* = hio/T is the dimensionless frequency. Note 
again that this spectrum is universal. The low-frequency 
cutoff scales with temperature in such a manner that 
hwp^-cci/T = Aa/J = const.. The value of the low- 
frequency cutoff is a fundamental constant given by Eq. 
(|6|) [correction due to /, Eq. ( |8a| ) is small and may be 
neglected] and independent of any parameters of the 
plasma. The universal, blackbody energy spectrum of 
radiation from the e"'"e~7 relativistic plasma is shown in 
Fig. n^. The Planckian spectrum is shown for comparison. 

At last we comment that the radiation from a relativis- 
tic fireball is always blackbody. Indeed, for a non-thermal 
radiation spectrum to occur, the linear size of the system 
should be i? < t / a-Tnpi{T ^ m^c^) < 10^^ cm (ctt is the 
Thompson cross-section and r < 1 is the optical depth 
of the plasma) which is microscopically small. 

Non-relativistic case — The universality of the pho- 
ton spectrum breaks down at plasma temperatures be- 
low T - 
(IPI) with 
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The dispersion relation is exactly Eq. 
replaced by its non-relativistic counter- 



To conclude, the thermodynamics of a gas of photons 
in the e^e~ plasma is studied. The deviation from the 
Planck distribution occurs because no photons with fre- 
quencies below the plasma frequency exist in the plasma. 
Equivalently, photons in the plasma acquire an effec- 
tive mass. In the ultra-relativistic plasma, T 3> ?TieC^, 
the effective photon mass scales as m^ — TA/c^ with 
A being the fundamental constant (independent of any 
plasma parameters) given by Eq. (^). Hence, the equilib- 
rium (blackbody) spectrum exhibits universal properties. 
The power-law scalings of thermodynamic parameters vs. 
temperature for "massive" photons are the same as those 
for the Planck distribution, but with different numerical 
pre- factors. In contrast, thermodynamics of photons in 
the non-relativistic plasma is not universal. 

The author is grateful to Ramesh Narayan and Abra- 
ham Loeb for useful and stimulating discussions. 
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where uq is the residual density of electrons at T 
(i.e., in absence of pair creation). 

A few estimates below show that the effect of the 
plasma frequency cutoff is usually small, unless a plasma 
is extremely dense and 'cold'. First, most of the explo- 
sion energy in cosmological, purely radiative fireballs (if 
any), — the 7-ray bursters, — is emitted when the sys- 
tem becomes optically thin for radiation. This occurs 
at r ~ 20 keV [0 when the e+e^-pair density drops to 
rigi ~ l/tTT^thin ~ 10^^ cm~^. The low- frequency cut- 
off is {hiOp/T) 
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10~ which is hardly observable. 
Second, distortions of the cosmic microwave background 
energy spectrum may be traced up to redshifts of z ~ 10^ 
pi. Taking the residual density of electrons no ~ lO^^ripi 
and temperature T r^ few keV, we estimate the cutoff 
{hiUp/T)(~.y^^ ~ 10~^. Third, in the center of the Sun 
uq ~ 10^^ cm"'^ and T ^ 1.4 keV which yields the cut- 
off (hup/T) 



Sun 



0.2; that is the cutoff is close to the 
Planckian thermal peak. The effect on the solar (and 
stellar) models will be addressed in future publications. 



where n^± = n^+ + n^- and contributions from heav- 
ier particles are neglected. When the temperature drops 
below TTieC^ the number density of e^e~-pairs decreases [1 
exponentially due to annihilation. Hence, we expect the 
non-relativistic cutoff frequency to be smaller than the 
ultra-relativistic. The density of charged particles is eas- 
ily calculated 0: 
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FIG. 1. The normalized blackbody equilibrium spectrum of 
radiation from the e^c~7 relativistic plasma and the Planck- 
ian spectrum vs. dimensionless frequency, cj, — hu/T. 



